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Numerical simulations of full QCD on anisotropic lattices provide a convenient way to study QCD thermo- 
dynamics with fixed physics scales and reduced lattice spacing errors. We report results from calculations with 
2-flavors of dynamical fermions where all bare parameters and hence the physics scales are kept constant while 
the temperature is changed in small steps by varying only the number of the time slices. The results from a series 
of zero-temperature scale setting simulations are used to determine the Karsch coefficients and the equation of 
state at finite temperatures. 
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1. INTRODUCTION 

The anisotropic formulation of lattice QCD has 
certain advantages when it comes to the study 
of the equation of state (EOS) at various finite 
temperatures. The finite temperature field the- 
ory has a natural asymmetry which makes the 
anisotropic approach useful to reduce the lattice 
spacing errors associated with the transfer matrix 
at less cost than is required for the full continuum 
limit Through the introduction of anisotropy 
on the lattice one can make the temporal lattice 
spacing, at, sufficiently small so that by varying 
only the number of time slices, TV*, the tempera- 
ture can be changed in small discrete steps. 

Our simulations of full QCD with two flavors 
of staggered fermions on anisotropic lattices are 
aimed at the study of the thermodynamic proper- 
ties of the quark-giuon system. We employ a fixed 
parameter scheme in which all the bare parame- 
ters of the simulation are kept constant and only 
the temperature is changed by varying iVt (from 4 
to 64). This approach separates temperature and 
lattice spacing effects and keeps the underlying 
physics scales fixed. 

The calculation of the EOS of the quark-gluon 
system involves derivatives of the bare param- 
eters with respect to the physical anisotropy 
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^ = as I at and the spatial lattice spacing Og. 
These "Karsch" coefficients are determined as 
a by-product of the zero-temperature calculations 
needed to choose the bare parameters. Once de- 
termined, the Karsch coefficients can be used for 
all temperatures since they depend only on the in- 
trinsic lattice parameters and not on Nt ■ This al- 
lows a straight-forward determination of the tem- 
perature dependence of the energy and pressure, 
again at fixed lattice spacing. With two or more 
slightly different values for at, a high-resolution 
sampling of temperatures can be investigated. 

2. THE ANISOTROPIC STAGGERED 
ACTION AND EOS DERIVATION 

Our calculations are based on the QCD action 
S^^S%^S%, where the gauge action is: 



^ P..,(a;)-fC2^F,t(a;) 



and the fermion action is: 
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To derive the EOS we start from the thermo- 



dynamics identities e(T) 
p(T) = T 
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Using the explicit form for Z and changing the 
independent variables with respect to which we 
differentiate to ^ and Os we get: 
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3. SIMULATIONS 

In our simulations we use the R-algorithm|^ 
implemented for two flavors of dynamical 
fermions. 

The zero-temperature runs from Table 1 are 
used to calculate the derivatives of the bare pa- 
rameters with respect to ^ and (the Karsch 
coefficients) involved in the EOS. 

The finite temperature runs in Table 2 and 3 
represent two separate sweeps through the tran- 
sition region for each of which the fixed bare pa- 
rameter scheme described above is applied. Fig- 
ure 1 shows the variation of (ipi/^) through the 
phase transition as we gradually change the tem- 
perature by varying only Nt for ^ = 4.0(1) and 
4.8(3). 
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16^x24x64 


1300 


5.695 


2.5 


0.025 
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16^x24x64 


1400 


5.725 


3.44 


0.025 
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16^x24x64 


3400 


5.6 


3.75 


0.025 
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16^x24x64 


3200 


5.3 


3.0 


0.008 


7 


16^x24x64 


3000 


5.29 


3.4 


0.0065 


8 


16^x24x64 


4300 


5.286 


3.427 


0.00394 



Table 1 

Parameters of zero-temperature calculations. All 
runs have dynamical vt = 1.0 except run 3 which 
has T^t ~ 1.2. All runs except 7 and 8 have valence 
i/f's of 0.8, 1.0 and 1.2. 
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3.0 
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5.3 


3.0 


0.008 
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Parameters of finite temperat 


are calculations 


c = 


4.0(1). All 


runs have dynamical vt = 1.0. 
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5.29 
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0.0065 
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10300 


5.29 


3.4 


0.0065 
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16^x12 


8800 


5.29 


3.4 


0.0065 
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16^x8 


3600 


5.29 


3.4 


0.0065 



Table 3 

Parameters of finite temperature calculations with 
^ — 4.8(3). All runs have dynamical ut = 1.0. 



4. KARSCH COEFFICIENTS 

We consider the bare parameters £,o, P, 
ruf and vt functions of the physical observ- 



ables Us, Rt = (m^/m^)*"°P°'-^' and 
= (TO2/m2)"P^"'^7(m2/m2)*™P°'-'^i and ex- 
pand those functions in Taylor series around the 
values of the bare parameters of a selected zero- 
temperature run. We fit the zero-temperature 
data to the linear parts of the Taylor expan- 
sions, where the fitting parameters then include 
the Karsch coefficients. The minimization in 
our case is implemented with the addition that 
we have to start with a guess for the standard 
deviations on the bare parameters and iteratively 
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improve this guess, yielding increasingly accurate 
values for the Karsch coefficients. Table 4 and 5 
show the numerical results for the Karsch coeffi- 
cients for expansion around the bare parameters 
of zero-temperature runs 6 and 7. 
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Table 4 










Karsch coefficients from 


expansion 


around run 6, Ta- 


ble 1. xl/D 


=(1.7, 1.0 , 


1.8, 0.7). 
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0.59(6) 


9.2(6.2) 
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0.7(1.8) 
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0.02(8) 


das 


-0.0050(8) 
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0.05(1) 


0.02(2) 
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0.06(3) 


-5.0(4.1) 


-0.8(7) 


0.8(1.2) 
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Table 5 
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=(1.4, 0.7 , 2.3, 0.9). 







5. EOS RESULTS 

Figure 2 shows the final result for the EOS. 
The errors on the pressure are significantly big- 
ger than the errors on the energy due to the large 
errors on those of the Karsch coefficients which 
are derivatives with respect to as. The compari- 
son with the free lattice theory (squares) gives an 
explanation of the prominent drop off of e and p 
in the high temperature sector — simply a con- 
sequence of the lattice high momentum cut-off. 
Our results are consistent with previous isotropic 
results 0. 

6. CONCLUSIONS 

We have studied the QCD thermodynamics 
using staggered fermions on anisotropic lattices. 
The fixed bare parameter scheme allows us to ex- 
plore the temperature dependence of energy and 
pressure with fixed physics scales. While this 
approach naturally reduces finite lattice spacing 
errors associated with at, the fixed lattice cut- 
off becomes important at increasing temperature. 
Including improvements to the spatial parts of the 
staggered fermion action would be a natural step 



to reduce the lattice artifacts for high tempera- 
tures. 
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Figure 1. The temperature dependence of (iljip) 

in the region of Tc- 



o E/e^g, ^ = 4.8(3) 

o p/PsB, ^ = 4.8(3)- 

♦ p/PsB, ^ = 4.0(1) - 

. e/Ejg, ^ = 4.0(1) _ 
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Figure 2. Energy and pressure in units of the 
Stefan-Boltzmann limit and a comparison with 
the free lattice theory (squares). 



